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ABSTRACT 

We have conducted a deep multi-wavelength analysis (0.15-160 /im) to study the outer region of 
the nearby galaxy M94. We show that the non-optical data support the idea that the outskirts of 
this galaxy is not formed by a closed stellar ring (as traditionally claimed in the literature) but by a 
spiral arm structure. In this sense, M94 is a good example of a Type III (anti-truncated) disk galaxy 
having a very bright outer disk. The outer disk of this galaxy contains ^23% of the total stellar mass 
budget of the galaxy and contributes ~10% of the new stars created showing that this region of the 
galaxy is active. In fact, the specific star formation rate of the outer disk (~0.012 Gyr -1 ) is a factor 
of ~2 larger (i.e. the star formation is more efficient per unit stellar mass) than in the inner disk. We 
have explored different scenarios to explain the enhanced star formation in the outer disk. We find 
that the inner disk (if considered as an oval distortion) can dynamically create a spiral arm structure 
in the outer disk which triggers the observed relatively high star formation rate as well as an inner 
ring similar to what is found in this galaxy. 

Subject headings: Galaxies: Evolution, Galaxies: Individual: Messier Number: M94, Galaxies: Indi- 
vidual: NGC Number: NGC 4736, Galaxies: Kinematics and Dynamics, Galaxies: 
Photometry, Galaxies: Spiral 



1. INTRODUCTION 

Recent years have seen a revolution in our understand- 
ing of the outskirts of spiral galaxies. On the one hand, 
UV emission in the disk of spiral galaxies has been discov- 
ered at very large radii (the so-called XUV disks Thilker 
et al. 2005; Gil de Paz et al. 2005; see also Donas et 
al. 1981), well beyond the optical radius. On the other 
hand, new studies have challenged the traditional picture 
of disks of spirals following an exponential radial surface 
brightness profile (Freeman 1970) down to a large radius 
where a truncation happens (van der Kruit 1979; van der 
Kruit & Searle 1981a, b). In fact, recent systematic obser- 
vations show an unexpected diversity of behaviors in the 
faintest regions of the surface brightness profiles. Three 
main categories or Types are now established: Type I, 
in which the disk does in fact show a simple exponen- 
tial form; Type II where the inner part of the disk shows 
an exponential fall-off followed by a steeper profile, and 
Type III, where the inner exponential profile is followed 
by a shallower, often exponential, profile (Pohlen et al. 
2002; Bland-Hawthorn et al. 2005; Erwin et al. 2005, 
2008; Hunter & Elmegreen 2006; Pohlen & Trujillo 2006). 
There is evidence (Pohlen & Trujillo 2006; Erwin et al. 
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2007) that these three types are related with the global 
morphology of the galaxies (i.e. Type II are more com- 
mon in later types), implying a suggestive link between 
the evolutionary patterns followed by the galaxies and 
the distribution of stars in their outskirts. 

Previous observational results suggest that the outer 
disks provide a direct view of disk assembly. For this rea- 
son, several theories have tried to explain the above vari- 
ety of behaviors in the outer regions of galaxies. Among 
them we can mention ideas from Yoshii & Sommer- 
Larsen (1989) who found that to explain Type I pro- 
files the time scales of both viscosity and star formation 
through the disks should be comparable (see also, Fer- 
guson & Clarke 2001 and Slyz et al. 2002). Also it is 
possible to find ideas based on star formation thresholds 
(Kennicutt 1989; Elmegreen & Parravano 1994; Schaye 
2004; Elmegreen & Hunter 2006), the maximum angu- 
lar momentum of the gas that is collapsing to form the 
disk (van der Kruit 1987; van den Bosch 2001), or sec- 
ular evolution driven by bars/spiral arms or clumps dis- 
ruption (Debattista et al. 2006; Bournaud et al. 2007; 
Roskar et al. 2008, Foyle et al. 2008) to explain the 
Type II morphology. To justify Type III, tidal stripping 
caused by minor galaxy mergers (Penarrubia et al. 2006; 
Younger et al. 2007), dark matter subhalos bombarding 
(Kazantzidis et al. 2008) or high eccentricity flyby by 
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a satellite galaxy (Younger et al. 2008) have been pro- 
posed. Although these theories qualitative explain the 
observed behaviors, a quantitative, detailed analysis of 
the observations that allow us to judge between the dif- 
ferent scenarios is missing. 

In this paper, we will focus on the nearby galaxy 
M94 (NGC4736). As we will show through the paper, 
this galaxy is an example of the brightest anti-truncated 
(Type III) outer disks that can be found in nearby galax- 
ies. The vicinity of the object (4.66 Mpc; Karachent- 
sev 2005) makes it a perfect candidate to explore in 
great detail the properties of its outer disk, and con- 
sequently, shed some light on the nature of the extended 
disk phenomenon, particularly, for the most brightest 
(anti-truncated) outer disk cases. 

To achieve this goal we have compiled data from UV 
(GALEX), optical (SDSS), near-infrared (2MASS) and 
infrared (Spitzer) to obtain the surface brightness distri- 
butions down to faint surface brightness in 17 different 
filters. The remainder of the paper is organized as fol- 
lows. In Section 2, we give a brief summary of the struc- 
tural characteristics of M94. Section 3 contains a detailed 
description of the data that we have used to conduct our 
analysis. The dependence with the wavelength of the 
outer galaxy structure is shown in Section 4. Section 5 
shows the surface brightness profiles of the galaxies in 17 
different filters from FUV to MIPS 160^m. The results 
of our analysis are explained in Section 6. Different for- 
mation scenarios are proposed to explain our results in 
Section 7. Finally, in Section 8 we summarize our main 
results. 

2. M94 (NGC4736) IN A NUTSHELL 

M94 is the closest early-type spiral galaxy with low 
inclination (see Fig. [1]). It is located in the Canes Ve- 
natici I cloud and was classified as (R)SA(r)ab-type ac- 
cording to de Vaucouleurs et al. (1991) and reclassified 
as (R)SAB(rs)ab-type by Buta et al. (2007). Its total 
apparent magnitude in the B-band is 8.54 mag, its R25 
is ~337" and it is claimed to be at an inclination of 35 
degrees (Erwin 2004). 

Up to five main regions have been identified in this sys- 
tem according to Bosma et al. (1977): a bulge R<15" 
(see also, Beckman et al. 1991), an inner ring of ongoing 
starburst activity (R«45"; van der Kruit 1976), an oval 
(disk/bar) stellar distribution down to intermediate dis- 
tance (R«220"; Beckman et al. 1991; Shaw et al. 1993; 
Mollcnhoff et al. 1995; Mulder 1995), a zone of low sur- 
face brightness (a gap), and what has been traditionally 
considered as an outer stellar ring around R«350" (see 
also, Buta 1988; Munoz-Tunon et al. 1989). Fig. [Uillus- 
trates the different regions of this galaxy as viewed by a 
very deep exposure obtained by our group (see details in 
the next section). 

M94 is also known for hosting the closest example of a 
low ionization nuclear emission line (LINER) 2 nucleus 
(see e.g. Roberts et al. 2001), although the origin of this 
emission as powered by a low-luminosity active galac- 
tic nucleus (LLAGN) is uncertain. Several observations 
supporting the hypothesis of a non-stellar origin for this 
emission are the detection of a strong non-thermal radio 
continuum source at the position of the nucleus (Turner 
& Ho 1994), the detection in UV using the HST of two 
bright point sources in the nuclear region with possi- 



ble bow shocks resulting from an hypothetical merger of 
two supermassive blackholes (Maoz et al. 1995), and the 
detection of a X-ray spectrum dominated by a Seyfert- 
like power-law continuum (Roberts, Warwick & Ohashi 
1999). On the other hand, favoring the stellar emis- 
sion hypothesis to explain the LINER emission are the 
large [Fell] /Pa/3 ratio (Larkin et al. 1998), and the Ha 
(Taniguchi et al. 1996) and FIR luminosities (Smith et 
al. 1991; 1994). 

3. OBSERVATIONS 

The analysis conducted in this paper is based on a 
multi-wavelength coverage of the galaxy. Consequently, 
the data comes from a variety of facilities and has been 
mostly obtained from archives. Below follows a brief de- 
scription of the data collected. 

3.1. Deep optical luminance imaging 

In order to explore in detail the structural character- 
istics of the outer region of M94 in the optical regime, 
we decided to take a very deep exposure of this galaxy 
(see Fig. Q]) using the Ritchey-Chretien 0.5-meter tele- 
scope of the BlackBird Remote Observatory (BBRO) sit- 
uated in the Sacramento Mountains (New Mexico, USA). 
The camera used was a Santa Barbara Instrument Group 
(SBIG) STL- 11 100 CCD camera, which has a FOV of 
27.7'x 18.2'and a pixel size of 0.45". The angular resolu- 
tion of this image is ^3.4". 

To maximize the photon collection our imaging con- 
sists on multiple deep exposures through a very broad 
non-infrared clear luminance filter (3500<A<8500 A), 
and a set of red, green and blue filters from the SGBI cus- 
tom scientific filters set 2 . The total luminance exposure 
time allocated to this image was 555 min. In addition, 
to decrease the effect of noise, and consequently increase 
the detection and enhancing of the faint structures we 
have applied a Gaussian blur filter (Davies 1990; Haral- 
ick & Shapiro 1992). The sigma of the gaussian kernel 
used on convolving the images is 1 pixel. This kernel was 
applied only to the outer regions to reveal more clearly 
this structure. As shown in a previous work (Martinez- 
Delgado et al. 2008) the faintest structures reachable 
with the BBRO telescope using the luminance filter with 
similar exposure times is >27.5 mag/arcsec 2 in the R- 
band. 

Since we lack a calibration of the optical luminance 
filter, our image is used only to show the faint structure 
of the outer disk at these wavelengths. 

3.2. GALEX UV 

The ultraviolet observations are obtained from the 
GALEX (Galaxy Evolution Explorer) mission (Martin 
et al. 2005). In particular, these data were retrieved 
from the Nearby Galaxy Survey mission (Bianchi et al. 
2003a,b; Gil de Paz et al. 2004; 2007). GALEX imag- 
ing benefits from a large field-of-view (1.25° diameter), a 
very low sky background and high sensitivity. We have 
used its two ultraviolet filters: a far-ultraviolet (FUV) 
filter centered at 1528 A(FWHM 269 A) and a near- 
ultraviolet (NUV) filter centered at 2271 A(FWHM 616 
A). The angular resolution of these imaging is ~4.6", and 

2 http://www.sbig.com/products/ 
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pixel size is 1.5". Absolute calibration uncertainties are 
~15% in both the far- and near- ultraviolet (Dale et al. 
2007). An image of M94 (NGC4736) in the ultraviolet 
filters can be seen in the left panel of Fig. [2] 

3.3. SDSS optical 

Optical data in 5 bands (u',g',r',i' and z') were ob- 
tained from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000) archive. SDSS images have a pixel size 
of 0.396"and typical angular resolution of Due to 

the large apparent size of the galaxy we had to create a 
mosaic of 8 individual SDSS frames (13.5'x9.8'each) to 
cover the total galaxy angular extension. An image of 
the galaxy in the SDSS optical bands can be seen in the 
central panel of Fig. O 

3.4. 2MASS near-infrared 

3 S , H and K s images were extracted from the 2MASS 
NASA/IPAC Infrared Science Archive (IRS A). This 
galaxy is included in the 2MASS Large Galaxy Atlas 
(Jarrett et al. 2003). This atlas already provides a sin- 
gle image (20' x 23.3') covering this galaxy composed of 
the individual 2MASS fields. The angular resolution is 
^3"and the pixel size is 1". 

3.5. Spitzer infrared 

IRAC (3.6/xm, 4.5/mr, 5.8/^m and 8fim) bands and 
MIPS (24/im, 70/im, and 160/zm) bands were obtained 
from the Spitzer Legacy program SINGS (Kennicutt 
et al. 2003). The pixel size is 0.75"for the IRAC 
(24.6'xl7.9') mosaics, and 1.5", 4.5"and 9"for MIPS 
24/tm, 70/j,m and 160/im mosaics respectively. Calibra- 
tion uncertainties are ~10% for IRAC data, and 4% 
(24/zm), 7% (70/im), 12% (160/im) for MIPS data (Dale 
et al. 2007). An image of M94 (NGC4736) in the Spitzer 
infrared filters can be seen in the right panel of Fig. [2l 
The angular resolution of these images is summarized in 
Table 1. 

4. OUTER GALAXY STRUCTURE WAVELENGTH 
DEPENDENCE 

The outer region of M94 (NGC4736) is claimed to be 
a stellar ring, as mentioned in Section 2. This is in fact 
the appearance of this structure when observed at the 
optical wavelengths (see the optical luminance image in 
Fig. [1] and the optical SDSS image in the central panel 
in Fig. [2]). The nature of this region, however, seems to 
be radically different when it is probed both in the UV 
and in the infrared (see left and right panels in Fig. [2]). 
In fact, at these wavelengths, the outer region is better 
described as an outer disk having spiral arms. 

Further support for the outer region being a disk with 
a spiral arm structure and not a closed stellar ring can be 
found from the deep HI imaging done by Braun (1995; 
his Fig. 8) and de Blok et al. (2008; their Fig. 80). In 
that Figures, a bright spiral arm is seen in the outer disk. 
This arm, that is not visible in the optical, corresponds 
to the brightest feature observed in the outer region on 
the Spitzer 24/im image (see Fig. ^ and it is also very 
clear in the UV imaging. Consequently, the non-optical 
data strongly points out in the direction of that the outer 
region of M94 is actually a disk with a spiral arm struc- 
ture and not a ring as has been traditionally claimed in 



the literature. On what follows we explore the properties 
of this outer disk in comparison to the different regions 
of the galaxy. 

5. SURFACE BRIGHTNESS PROFILES 

In order to derive surface brightness profiles in a con- 
sistent way we use the same set of ellipticity and posi- 
tion angle in all our bands. Following Pohlen & Trujillo 
(2006) we explore which value of the ellipticity and po- 
sition angle (PA) are a good representation of the outer 
disk. To remove as much as possible potential effects re- 
lated to star formation and dust extinction we use the 
3.6/im image as a proxy. Briefly, the method works as 
follows, we use the IRAF task ellipse (STSDAS package) 
on the sky substracted and masked (from all non-galaxy 
component) image. The center of the ellipses is fixed by 
a gaussian fit to the bright nucleus of the galaxy. We 
used a logarithmic radial sampling with steps of 0.03 to 
increase the S/N especially in the outer parts. Iterative 
3cr rejection along the ellipse is applied to minimise the 
influence of cosmic rays or remaining foreground stars. 
The free ellipse fit (fixed center, free ellipticity and posi- 
tion angle) is used to determine the best set of ellipticity 
and PA describing the outer disk. These values are taken 
at the radius where the mean flux of the best fitted free 
ellipse reaches the value of the standard deviation of our 
background sky. This limit ensures enough S/N to fit a 
free ellipse but is small enough to be in the radial region 
dominated by the outer disk. 

Using 3.6/^m image, we find that the outer isophotes 
are well described by an ellipticity of 0.14 and a PA of 
122.2 degrees (anticlockwise from North). This ellipticity 
implies an inclination of ^8 degrees for the galaxy, far 
away from previous reported values of 35 degrees based 
on the shape of the inner (oval distortion) disk. The PA 
is also different from the one reported previously of 108 
degrees (Sofue et al. 1999) for the same reason. 

The surface brightness profiles in all the filters were ob- 
tained using a task custom-made using IDL. Bright stars 
were identified and masked. We measured the median 
intensity and rms along elliptical apertures with fixed el- 
lipticity and position angle. Our code was tested against 
the IRAF task ellipse producing very similar results. 

A critical issue for obtaining accurate surface bright- 
ness profiles at the low surface brightness levels of the 
galaxies is the accuracy of the sky substraction. Al- 
though the archive data used in this work is sky sub- 
tracted we repeated the "ellipse method" explained in 
Pohlen & Trujillo (2006) to improve the sky determina- 
tion. In brief, the idea is to expand the surface brightness 
estimation using the ellipticity and PA obtained from the 
outer disk. By plotting the flux at these ellipses as a func- 
tion of radius it becomes clear at which radial distance 
the ellipses leave the galaxy and enter the background, 
flat flux (noise), region (see Fig. 2 of Pohlen & Trujillo 
2006). The mean value and standard deviation of the 
fluxes within the above radius and the radius extended 
by 20% are used as the final sky values. This sky value 
is subtracted to the original image and the process of 
estimating the surface brightness is repeated. Our sur- 
face brightness estimation (given in the AB system) are 
shown on Fig. [3l 

Each image has a critical surface brightness (/J. C rit) be- 
yond that the reported surface brightness profiles are not 
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reliable. Using SDSS imaging, Pohlen & Trujillo (2006) 
defined this [i c nt as the value where the profiles obtained 
by either over- or undersubtracting the sky by ±l<r-sky 
start to deviate by more than 0.2 mag. This value is 0.54 
mag above the limiting surface brightness (fiu m ) of the 
image calculated as the 3<r error of the sky value determi- 
nation. Table [1] compiles the set of fi C rit for the different 
bands used in this work. The surface brightness are cor- 
rected for Galactic extinction (Schlegel et al. 1998) using 
E(B-V)=0.018 mag. No attempt was made to correct the 
surface brightness measurements for internal extinction, 
since no unique recipe is available to do this. Conse- 
quently, this will introduce an additional uncertainty on 
the surface brightness estimation. In addition, the galaxy 
studied here is a fairly face-on system (the measured el- 
lipticity implies an 8 degree inclination) so the expected 
correction is only very small. In fact, if we adopt a pop- 
ular parameterization of the effects of internal extinction 
(see e.g. Cho & Park 2009): A A =7.\log 10 (a/b) where 7a 
is a constant that depends on the observed wavelength 
but which value is ~1, and a/b is the axis ratio, then in 
our case, log 10 (1.16)=0.064. So the expected corrections 
are really minor. 

Following the estimations of fi C rit we find that, with 
the exception of the shallow 2MASS near-infrared im- 
ages (which do not allow us to study the outer disk), we 
can choose the following outest radius (430") as a com- 
promise among the different profiles depth. From now 
onwards, consequently, we will divide the galaxy into the 
following zones: bulge region (0-75"), inner disk region 
(75-200"), and outer disk region (200-430"). 

6. RESULTS 

6.1. Surface brightness profile shapes 

The observed profiles coupled with the galaxy mor- 
phology suggests classifying this galaxy as a Type III 
(Erwin et al. 2005). In fact, one of the most interesting 
features of the surface brightness profiles of M94 is the 
unchanging nature of the global structure of the galaxy 
from the UV to the far infrared. This implies that stars of 
all ages as well as dust are similarly distributed through- 
out the galaxy structure. The similarity between the UV 
and FIR profiles is more or less expected as both wave- 
length ranges are tracing the most active star formation 
regions of the galaxies. In relation to the optical and 
NIR data (although we warn that for these last data we 
lack the information in the outer region) which trace the 
intermediate-old stellar component, we think that the 
fact that they also share the shape of the young stellar 
component implies that the galaxy, globally, is actively 
forming new stars through all the disk and that this ac- 
tivity has been maintained similarly in the last Gyrs. 
As we will show later in the text, this scenario is sup- 
ported by the specific star formation rate distribution of 
the galaxy which is relatively constant within a factor 
of 3 through the whole galaxy and by the fact that the 
gas exponential infall time scales are similar along the 
galaxy. 

An exception is seen for the 70 and 160 /im MIPS 
bands, where we see the inner star forming ring is al- 
most invisible. However, that is clearly the effect of the 
larger PSF at these wavelengths (cf. Table 1). 

One of the most interesting feature of the surface 
brightness profile of M94 is the brightness (~23.5 



mag/arcsec 2 in r'-band) at which the outer disk appears. 
Normally, for Type III galaxies, the break occurs at 
faintest surface brightness. For example, Pohlen & Tru- 
jillo (2006) find a typical surface brightness of ~24.7±0.8 
mag/arcsec 2 in the r'-band, whereas Erwin et al. (2008) 
find ~24.2±0.5 mag/arcsec 2 in the R-band. This makes 
M94, in comparison, a Type III galaxy having an outer 
disk on the very bright side. However, M94 is not an 
unique example within the family of Type III objects, 
in fact, it is possible to find other galaxies with similar 
characteristics. For example, M77 (NGC1068) is strik- 
ingly similar to M94 (see its profile on Pohlen & Trujillo 
2006). In the case of M77, the outer disk starts at ~23 
mag/arcsec 2 , and morphologically also, both have simi- 
lar features. For this reason, we think that both of these 
two galaxies could be considered as examples of Type 
III galaxies having the brightest (<23.5 mag/arcsec in 
r'-band) anti-truncated outer disks. However, a note of 
caution should be stated about the Type of the galaxy 
M94. The Type III classification is based on the assump- 
tion (following the traditional picture of this galaxy) that 
the inner region (75-200") is a true disk. If this region 
is, however, considered as an oval distortion (i.e. a struc- 
ture with a bar-like nature and not a real disk) , then the 
genuine disk of the galaxy is the outer region (>200") 
and the galaxy should be better considered as a Type I. 

6.2. Stellar mass surface density profile 

At any given wavelength it can be shown that the sur- 
face brightness profiles /x A can be converted into stellar 
mass surface density profiles as follows: 

log(M/(M pc 2 )) = log(M/i) A -O.4( AiA -m Qbs . . A )+8.629 

(1) 

with (M/L)\ the stellar mass-to-light ratio in solar 
units and m a b s ,Q,\ the absolute magnitude of the Sun 
at a given wavelength A. Lacking a deep near-infrared 
image, we have chosen the SDSS r-band to estimate the 
M/L ratio. This band is selected because it is the reddest 
of the deepest SDSS bands. To estimate the (M/L) r ratio 
we have used the (g-r) color following the prescription 
given by Bell et al. (2003). We add 0.15 dex to the M/L 
obtained in Bell et al. (2003) to normalize the predicted 
M/L to a Salpeter (1955) IMF. The (g-r) color profile, the 
(M/L) r and the stellar mass density profile are shown in 
Fig. [H The minimum observed both in color and (M/L) r 
is associated with the nuclear stellar ring. The shape of 
the stellar mass density profile is very similar to the shape 
of the reddest (stellar origin) surface brightness profiles 
we have. The stellar mass density at which the outer 
disk starts is ~60 M Q /pc 2 . As expected due to M94's 
bright outer disk, this value is higher than the mean value 
measured for late-type Type III galaxies (~10 M Q /pc 2 ; 
Bakos et al. 2008). 

Using the stellar mass density profile based on colors, 
one can calculate a rough estimate of the total stellar 
mass of the galaxy: ~6.5 x 10 10 M Q . In each of the galaxy 
regions 3 the amount of stellar mass we measure is: bulge 

3 We would like to emphasize that the term regions used here 
(as well as in the rest of the paper) should not be confused with 
the actual components of the galaxy (i.e. bulge, inner disk, outer 
disk, etc). This distinction is relevant when estimating quantities 
like the amount of mass since the real components of the galaxy 
can indeed overlap. 
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region (3.2xl0 10 M ), inner disk region (1.7xl0 10 M Q ) 
and outer disk region (1.5xl0 10 Mq). Due to the signif- 
icant brightness of the outer disk, its contribution to the 
total stellar mass budget of the galaxy is quite significant 
(~23%). This is a factor of —2 larger than the average 
late-type Type III outer disk (Bakos et al. 2008). 

6.3. Star formation rates 

To assess the origin of the outer disk structure a key 
element is to measure its star formation rate. This value 
can give us clues on whether this region of the galaxy 
is still forming (i.e. creating new stars) or whether it is 
an old structure. Following Munoz-Mateos et al. (2007), 
the star formation rate (SFR) can be computed from the 
apparent magnitude (corrected for internal extinction) 
in the FUV (AB system) using the calibration given by 
Kennicutt (1998): 

log(£F J R)(Af yr -1 ) = 2 logd(pc)-0.4F[/U-9.216 (2) 

This can be transformed to a SFR per surface area as 
follows: 

logiSFR^iMayr-'pc- 2 ) = -0.4/x Fi7V + 1.413 (3) 

To estimate the internal extinction Apuv,in hi the 
FUV, we follow the empirical relation between the TIR- 
to-FUV ratio (i.e. the total luminosity in the infrared 
over the luminosity in the FUV), log(Fr/fl/ 'Fpuv) and 
the observed (FUV-NUV) color derived by Boissier et al. 
(2007). Once log{FTiR/ 'Fpuv) is obtained the attenu- 
ation profiles are computed using the fit of Buat et al. 
(2005): 

Apuvin = -0.0333V 3 + 0.3522V 2 + 1.1960V + 0.4967 

(4) 

with X=log(FTiR/FFuv)- The median value of 
A-Fuv,in within the inner 430" is 1.7 mag. This translates 
(using A K =0.0465A F;7 y; Munoz-Mateos et al. 2007) to 
an inner median extinction in the K-band profile of 0.079 
mag. 

Using Eq. [2] we can estimate the total SFR of the 
galaxies as well as the SFR at the different components. 
We find that the total SFR for M94 is 1.04 Mgj/r" 1 
which is distributed across the different regions of the 
galaxy as follows: bulge (0.75 M Q 2/r _1 ), inner disk (0.14 
M yr _1 ) and outer disk (0.15 M j/r _1 ). The amount of 
stars the outer disk is forming is comparable to the inner 
disk and it is a significant (~ 10-15%) fraction of the total 
star formation rate of the galaxy. This clearly shows, 
that at least for M94, the outer disk is a very active 
region of the galaxy. This activity is directly visible when 
looking at the UV and IR images shown in Fig. [2j where 
the unobscured young stars as well as the hot dust are 
clearly seen in form of spiral arms in the outer region of 
the galaxy. 

Another important parameter to measure is the spe- 
cific star formation rate (sSFR), i.e. the amount of star 
formation per unit stellar mass. This value informs about 
the efficiency of the star formation in the different regions 
of the galaxies. The sSFR can be derived from the com- 
bination of Eqs. [1] and [3] to give: 



logisSFfyiyr- 1 ) = -0A(FUV - m x ) - log(M/£ A ) 

-0.4™ abS!O , A - 7.21^5) 

Both the SFR and sSFR radial distribution are 
plotted in Fig. [5] Our estimates agree very well 
by those produced by Boissier et al. (2007) for 
the star formation rate density (their Fig. 9.39) 
of this galaxy: ~lO _8 M pcr 2 2/r -1 at 50"(bulge), 
— 5x 10~ 9 M Q pc~ 2 yr~ 1 at 150"(inner disk) and 
~4xlO _10 Af Q pc~ 2 yr" 1 at 300" (outer disk). Munoz- 
Mateos et al. (2007) have noted that nearby spirals show 
a slightly positive increase of the sSFR from inside-out. 
This result also applies in the particular case of M94, 
where we find that the outer disk has a larger sSFR 
(-0.012 Gyr" 1 ) than the inner disk (-0.006 Gyr^ 1 ). 
In other words, per unit stellar mass, the outer disk is 
more active (efficient) forming stars. 

6.4. Spectral Energy Distributions 

Using the fluxes obtained in the 17 filter bands we have 
created the Spectral Energy Distributions (SEDs) rang- 
ing from 0.15-160/im in the different regions of the galaxy 
(see Table [2]). The SEDs across the galaxy are shown in 
Fig. El 

To extract some relevant physical parameters from the 
SEDs we have used synthetic templates computed with 
the code GRASIL (Silva et al. 1998). This code couples 
the spectral (and chemical) evolution of stellar popula- 
tions with the radiative transfer through a dusty inter- 
stellar medium. Consequently, the models consistently 
account for the emission from stars as well as absorption 
and emission from dust in galaxies. GRASIL allows the 
variation of a large set of free parameters but here we 
concentrate only in the small subset that mostly affect 
the shape of the SED (Panuzzo et al. 2007). We describe 
briefly which parameters were left to change during the 
fitting of the SEDs (see Silva et al. 1998 for a detailed 
description of the code): 

1. Total infall mass, Mj„/. This is the amount of pri- 
mordial gas that has infalled to form the galaxy. 
We have allowed the following range of variation 
depending on the region of the galaxy which is fit- 
ted: bulge (0.8 to 5xl0 10 M Q ), inner disk (0.6 to 
3x 10 10 M o ), outer disk (0.6 to 3x 1O 1O M ) and the 
full galaxy (1 to 7xl0 10 M o ). 

2. Exponential infall time scale, 77,. Gas is assumed 
to be continuously infalling at a rate proportional 
to exp (— t/rb). n was allowed to change from 0.1 
to 10 Gyr. 

3. Efficiency of Schmidt-type law, v. The SFR is as- 
sumed in GRASIL to follow a Schmidt-type law 
of the form SFR(t)=z/M s (t) fe , where M s is the gas 
mass at any time, k the exponent of the Schmidt 
law which we fixed to 1 , and v is the free efficiency 
parameter varying from 10 -3 to 4. 

4. Fraction of molecular mass to total gas mass 
(M mo ;/M gas ). We allow this parameter to range 
from 0.1 to 0.99. 
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5. The escape timescale of newly born stars from their 
parent molecular clouds, t esc . We let this vary from 
1(T 3 to 1CT 1 Gyr. 

6. The mass of the molecular clouds, Mmc- The mass 
of the molecular clouds together with the radii of 
the molecular clouds is related to the optical depth 
of the molecular clouds, tmc °c MmcI t \ic- We 
fixed rMC=14 pc and let Mmc vary from 10 3 to 
10 7 M Q . 

To allow a comparison with other works the adopted 
IMF is a Salpeter IMF in the mass range from 0.1 to 100 
Mq . Due to the low inclination of our galaxy we use the 
output of the code with zero inclination. We assume that 
the galaxy has a present age of 13 Gyr. 

For the SED of each region we have computed the min- 
imum x 2 ■ Due to the large number of free parameters (6 
in our case) we have avoided a brute force approach that 
implies to compute a large set of templates. For exam- 
ple, assuming that ten values per parameter are desired, 
lQN par templates should be created. In current powerful 
desktops every GRASIL template takes ~30 seconds to 
be computed. Consequently, the amount of time required 
in a brute force approach is prohibitive. To overcome this 
difficulty, we have applied a Markov chain Monte Carlo 
technique based on the Metropolis algorithm (Metropo- 
lis et al. 1953; Neal 1993). The details of the code that 
we have used are in Asensio Ramos et al. (2007). 

The best fit parameters together with the 1 a error are 
shown in Table [31 Together with these values GRASIL 
also provides the present SFR and the optical depth of 
molecular cloud at 1 /im tmc and the amount of mass 
in gas and stars. 

We can conduct now a comparison between the 
stellar mass obtained using the (g-r) color and the 
ones obtained using the GRASIL code. Compar- 
ing the stellar mass determined with GRASIL with 
the ones derived in Section 6.2, following Bell et al. 
(2003), we find them systematically (~0.7) lower. 
However, it is worth stressing that within the error 
bars both estimations agree. From the above com- 
parison we think it is safe saying that the amount of 
stellar mass contained within the different regions of 
the galaxy is as follows: (2.8<M*<3.2) x 1O 1O M 
(bulge), (1.3<M*<1.7)xlO lo M (inner disk), 
(O.9<M*<1.5)xlO lo M (outer disk), and 
(4.2<M*<6.5)xlO lo Af (total). In relation to the 
SFR, we find that GRASIL and the prescription pre- 
sented in Munoz-Mateos et al. (2007) agree very well 
within the error bars. 

The GRASIL prediction in relation to the amount of 
present gas (i.e. not transformed into stars) in the differ- 
ent regions of the galaxy can be also compared with the 
recent measurements of HI (Walter et al. 2008). These 
authors find a total amount of gas of 4x 1O 8 M for M94. 
This number is in excellent agreement with the GRASIL 
prediction (4.5±l.O)xlO 8 M . 

From the GRASIL results one can extract the following 
information: 1. The efficiency of the star formation, as 
characterized by the v parameter is very similar through 
the whole structure of the galaxy, with no significant 
changes from the bulge to the outer disk. 2. The mass of 
the typical molecular clouds where the star formation is 



taking place seems, however, to be much smaller in the 
inner region of the galaxy than in the outer region. This 
could be related to the much more violent conditions 
(influence of the active nucleus) close to the center of 
the galaxy preventing the formation of massive molecular 
clouds. 3. The exponential infall time scale of the gas 
seems to vary also across the galaxy. The gas appears to 
take more time to be deposited in the bulge than what 
requires to be placed in the outer disk. This could be 
understood if the gas requires, as expected, a large loss 
of angular momentum to feed the bulge than the inner 
disk. 

7. ORIGIN OF THE OUTER DISK 

In the following we will address the origin of M94's 
bright outer star-forming disk. We will explore the two 
obvious possibilities of it having an external versus inter- 
nal origin. 

7.1. Merger origin 

It has been proposed (Penarrubia et al. 2006) that 
extended galactic disks can be form by a tidally disrupted 
dwarf galaxy in prograde orbit that is coplanar with the 
host galaxy disk. As a result of this merging, the stellar 
debris ends up resettling into an extended exponential 
disk in dynamical equilibrium. Depending on the mass of 
the satellite, orbit eccentricity and inclination, the stellar 
debris exhibit a rotation that could range from 30-50 
km/s lower than the host circular velocity up to a factor 
of ^2 smaller. 

Favoring a merger origin for the outer disk of M94 is 
the fact that its rotation curve is declining from the inner 
region (~3 kpc) to the outer region (~9 kpc) by ~30% 
(Mulder 1995; de Blok et al. 2008). However, a common 
concern with the merging scenario (particularly if the 
merger has a large mass ratio) is the possibility that the 
accretion destroy the inner disk (Hopkins et al. 2009). In 
our case, the outer disk of M94 contains a large fraction 
(^23%) of the total mass of the galaxy. If all the outer 
disk were entirely the result of a disrupted dwarf the 
implied merger ratio would be as high as 1:4. 

Another problem with the merger scenario, if we want 
to explain entirely the outer region of M94 as a tidally 
disrupted satellite, is the very different ellipticity the 
inner disk and the outer region has. In the merging 
hypothesis the ellipticities of both components can be 
interpreted as a reflect of the different inclination they 
have. In fact, transforming the observed ellipticities to 
inclinations, we find that the inner region will have an 
inclination of ^35 degrees whether the outer disk will 
have only ~8 degrees. Merging simulations of massive 
satellites with disk galaxies whereas the disk survived 
the encounter only found a maximum tilt of ~7 degrees 
for the outer regions (Velazquez & White 1999). This is 
much smaller than what is observed for this galaxy. 

Even in the case that the origin of the outer region of 
M94 could not be explained by a major merger it could 
be well the case that an interaction of M94 with another 
galaxy could be the responsible of triggering both the 
star formation we see in this part of the galaxy as well as 
the many morphological structures this galaxy has: inner 
ring, inner disk/oval distortion, gap, outer spiral arms, 
etc. The interaction possibility is however disfavored by 
the observations. Although M94 is the main member of 
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the Canes Venatici I Group, it seems unlikely that it has 
suffered any interaction with another member of that 
group recently. In fact, such group has been found to 
be only weakly gravitationally bound. Moreover, most 
of the galaxies in this group appear to be moving with 
the expansion of the universe, which will prevent nearby 
(repeated) encounters affecting M94 (Karachentsev et al. 
2003). 

7.2. Secular evolution 

We also explore whether the plethora of structures we 
see in M94 as well as the enhanced star formation in its 
periphery could be created assuming a perturbation of 
internal (secular) origin. It is well known that bars can 
be produced via dynamical instabilities in galactic disks. 
A particular form of bars is what is known as oval dis- 
tortion which are weak bars also called "fat bars" that 
often resemble normal slightly inclined disk having even 
spiral arms (see e.g. Jogee et al. 2002). As mentioned in 
Sec. 2, the inner disk of M94 has been sometimes con- 
sidered as an oval distortion more than a truly inclined 
disk. This seems to be a reasonable hypothesis since the 
inclination of the galaxy, as derived in the UV and IR 
imaging for the outer disk, is lower than what it would 
be derived from the shape of the inner disk. It is worth 
stressing that this conclusion suggesting that the inner 
region of M94 is an oval distortion more than an truly in- 
clined disk can not firmly established on the basis of the 
optical images alone (which were the main source of in- 
formation used in previous works), but, both the UV and 
IR images presented in this paper clearly show that there 
is a continuation in the spiral arms from the inner (oval) 
region to the outer disk. If both structures were decou- 
pled, as suggested by the apparent different inclination 
and explained in the merger scenario, this prolongation 
of the spiral arms would be extremely unlikely. Under 
the hypothesis that the inner region is an oval distortion 
the problem with the different inclinations among the 
inner and the outer disk of M94 is solved since the oval 
distortion is not inclined but has a different cllipticity 
than the outer disk because is intrinsically oval and not 
circular. 

If the inner region is an oval distortion, this structure 
could be a serious candidate to explain the spiral arm 
structure we see in the outer disk (Binney & Tremaine 
1987, and references therein), and consequently to be 
the ultimate responsible of the enhanced star formation 
we measure in this part of the galaxy. Oval distortions 
are known to create spiral arms that can be sustained 
for long time in isolated galaxies (e.g., Lindblad 1960; 
Toomre 1969; Sanders & Huntley 1976; Athanassoula 
1980). In fact, oval distortions have been considered as 
important as bars when considering the secular evolution 
of a disk galaxy (Kormendy & Kennicutt 2004). 

We have created a toy model to explore whether the 
observed oval distortion can produce spiral arms similar 
to what we see in the outer disk. Our simple model is 
composed of a bulge, an oval distortion, an outer disk 
and a dark matter halo. To characterize the different 
components of the galaxy we use different analytical for- 
mulas that we fit to the stellar mass density profile in 
order to have realistic parameters (see Fig. [7]). Consid- 
ering the oval distortion as a very mild bar, we use the 
analytical formula by Ferrers (1877) to describe it: 
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where m 2 = f^- + p- + fr , and a is the major axis, b the 
minor, and c the vertical axis, taken here as 500 pc, with 
a > b > c. This bar model is very popular in the litera- 
ture (see e.g. Sanders & Tubbs 1980; Papayannopoulos 
& Petrou 1983; Schwarz 1985; Athanassoula 1992; Patsis 
2005, etc). From the fit to the stellar surface mass den- 
sity we obtain: a=4.9 kpc, b=4.1 kpc and a total mass 
for the oval distortion equal to M a d = 1.3 x 1O 1O M0. 

The remaining components of the galaxy (the bulge, 
the outer disk and the dark matter halo) have been mod- 
eled using the Miyamoto & Nagai (1975) analytical mod- 
els. The Miyamoto-Nagai potential is described by the 
equation: 



®m(R, z) 



GM 



^R 2 + (a + Vz 2 + b 2 ) 2 
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The parameters retrieved from the fit to the stel- 
lar mass surface density are: bulge (M=2.0x 10 10 M Q , 
a = 0.05 kpc and b — 0.6 kpc); outer disk 
(M=2.8xlO lo M , a = 6.5 kpc and b = 0.5 kpc) and 
halo (M=2O.OxlO lo M , a = 0.0 kpc and b = 30 kpc). 
Note that the outer disk in this model has a factor of 
~2 more stellar mass than what we have measured be- 
fore using the outer region (200-430") of the stellar mass 
profile. This is explained because in the model the stel- 
lar mass density profile is extrapolated in- and outwards, 
to zero and infinity respectively. Once the galactic po- 
tential is set, we can start running a simulation and see 
how a gaseous disk responds to the potential. Before do- 
ing this we can use the linear approximation to calculate 
the angular frequency, and therefore estimate the radii 
of the main resonances (Fig. [8]) associated to the grav- 
itational potential we have created. This is important 
for this galaxy since several resonant structures like the 
inner ring are observed. 

One of the key issues when considering bars, or oval dis- 
tortions, and their corresponding resonances is the pat- 
tern speed, i.e., how fast the bar rotates as a solid body. 
The pattern speed determines, many different structures, 
such as the radii of the different resonant rings, or their 
actual existence. In the case of our study, the variation 
of the pattern speed, also determines the appearance of 
the spiral arms and their shape. A precise study of the 
spiral arms and their origin (Romero-Gomez et al. 2007), 
establishes that when increasing the pattern speed in a 
given barred model, the spiral arms become more tightly 
wounded. 

Several pattern speeds have been proposed to 
explain different features in M94. We started 
with O p =56 km s -1 kpc -1 proposed by Mulder & 
Combes (1997) and decreased this value down to 
fi p =30 km s -1 kpc" 1 . We found that a model that is able 
to reproduce the shape of the outer spiral arms as well 
as the inner ring of the galaxy has f2 p =38 km s -1 kpc -1 . 
This number agrees nicely with previous studies, such as 
Waller et al. (2001), where they studied the central part 
of the galaxy and obtained f2 p =35 km s -1 kpc -1 . With 
the pattern speed that we are using, Rch = 5.2 kpc, and 
consequently R C 7fl/R( )ar =1.04. In that sense, the oval 
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distortion that we have created would be considered as 
a fast bar (Aguerri et al. 2003). In addition, our pattern 
speed produces an inner ring at a radius of ~1 kpc, at 
the same distance than the one observed. 

Once the gravitational potential and the pattern speed 
are fixed we run our simulation to explore the response 
on a pure gas disk. We use the code FTM 4.4 (updated 
version) from Heller & Shlosman(1994), which allow us 
to study the dynamical response of the gas based on 
SPH. We use 150000 collisional particles to simulate the 
gas, isothermal and non self-gravitating. Initially the 
gaseous particles are distributed following a Miyamoto- 
Nagai density distribution. In Fig. [9l left panel, we show 
a snapshot of our simulation showing clearly the spiral 
arms, as a result of the oval distortion. In addition, the 
inner ring that is observed in the real galaxy can be 
seen. To make a more quantitative comparison of our 
toy model with the observations, in Fig. [9l right panel, 
we compare the density profile of the gas in our simula- 
tion with the FUV surface brightness profile. The peaks 
on the FUV profile represents the places where the star 
formation is enhanced. We can see that there is a rea- 
sonable agreement between the peaks in the FUV and 
the position of the peaks where the density of the gas in 
the model is higher. 

According to our simulation, if the oval distortion sur- 
vives for long time (Athanassoula & Misiriotis 2002, 
Martinez- Valpuesta, Shlosman & Heller 2006) the spi- 
ral arm and the ring will also survive. However, a slow 
secular evolution is expected in the radius of the oval 
structure and the pattern speed. In addition, we should 
take into account the slowly consumption of the gas by 
the star formation (which could be somehow counterbal- 
ance by the new gas infall), that will also decrease the 
amount of gas in the spiral arms and oval distortion. 

In a following paper (Martinez- Valpuesta et al., in 
preparation) we plan to increase the resolution of our 
hydrodynamical simulation to explore in more detail the 
inner region of the galaxy. In particular, we will explore 
the velocity shear field of the galaxy to probe whether 
the dynamics linked to the oval distortion can help to ex- 
plain the GRASIL findings in relation to the efficiency, 
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molecular cloud mass and infall time scale found from 
the stellar population analysis. Our aim is to find a con- 
sistent model able to reproduce at the same time the new 
HI kinematic data (Blok et al. 2008) in the outer region. 

8. CONCLUSIONS 

M94 is the closest example of an early-type spiral 
galaxy with low inclination. It has been traditionally 
considered (due to its optical appearance) as a galaxy 
with a bright outer ring surrounding its disk. The mul- 
tiwavelength analysis that we have conducted in this pa- 
per show that this view is incorrect. Observed in the UV 
and in the IR the galaxy reveals that its outer region has 
a disk with a spiral arms structure nature. This outer 
disk is very active forming stars and contains ~23% of 
the total mass of the galaxy. We have explored what 
is enhancing this star formation activity in the galaxy 
periphery and find that a likely candidate is the dynam- 
ical effect of the inner oval distortion on the rest of the 
galaxy. Although an external origin of the extended disk 
(e.g. due to the accretion of a satellite galaxy) can not 
be ruled out with the existing data, the oval distortion 
hypothesis explains more naturally the existence of an 
inner ring and the enhanced star formation activity in 
the outer disk of M94. In addition, it is likely that the 
flux of gas that the oval distortion produces towards the 
center of the galaxy would be related to the active nature 
of the bulge of this galaxy. 
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Fig. 1. — Luminance clear-filter very deep image of M94 (NGC4736) obtained with the BBRO 0.5-meter telescope. The total exposure 
time of the image was 555 min. The image size is 27.7'xl8.2'. North is up and East is left. For illustrative purpose, a color image of the 
central part of the galaxy obtained with the same telescope has been superimposed on the inner (saturated) disk region of the galaxy. It 
is worth stressing that the traditional ("textbook") picture of M94 is just what it is shown in color here. 




Fig. 2. — Panel showing M94 (NGC4736) galaxy through the ultraviolet, optical and infrared wavelengths. Left panel shows M94 as 
viewed by the GALEX bands (blue: FUV; yellow: NUV). Middle panel shows M94 as viewed by the SDSS bands (blue: g-band; green: 
r-band; red: i-band). Right panel shows M94 as viewed by the Spitzer bands (blue: 3.6/Ltm; green: 8fim; red: 24/xm). In all the cases North 
is up and East is left. The area shown (18'xl5') is selected for being the common area covered by the set of facilities used in this work. 



Panchromatic view of M94 (NGC4736) 
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Fig. 3.— Surface brightness profiles (AB system) of M94 (NGC4736) in 17 bands ranging from the FUV to MIPS 160^tm. 



TABLE 1 

Data quality and galactic extinction correction 



Band ^ crit (mag/arcsec 2 ) R crlt (") PSF FHWM (") A A (mag) 



FUV 


31.5 


550 


1.5 


0.134 


NUV 


31.3 


650 


4.7 


0.119 


u 


28.2 


660 


1.5 


0.093 


g 


28.2 


710 


1.2 


0.070 


r 


27.5 


730 


1.0 


0.050 


i 


27.2 


690 


1.1 


0.037 


z 


26.7 


670 


1.0 


0.030 


.1 


24.1 


205 


3.0 


0.017 


H 


23.8 


210 


2.9 


0.010 


K 


24.3 


210 


3.0 


0.007 


3.6/jm 


26.6 


660 


1.7 


0.004 


4.5/^m 


25.8 


630 


1.9 


0.001 


5.8/im 


25.0 


360 


2.0 




8.0/mi 


26.0 


415 


2.2 




24/^m 


25.5 


360 


5.5 a 




70^m 


23.2 


410 


15.8 a 




160/xm 


21.5 


450 


34.6 a 




Note. — 


Last column, Aa, 


is the magnitude 


correction applied to corre 



for Galactic extinction (Schlegel et al. 1998) using E(B-V)— 0.018 mag. 
a Values obtained from a Gaussian fit to theoretical PSF models constructed 
by the Tiny Tim Point Spread Function modeling program (developed for the 
Spitzer Science Center by John Krist; STScI) to a blackbody source of 50K 
temperature ( http: / / ssc.spitzer.caltcch.edu/mips/psf.html) . 
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Fig. 4. — Left panel: (g-r) color and (M/L) r profiles of M94. Right Panel. Stellar mass density profile of M94 (NGC4736) obtained using 
the r-band profile and (M/L) r from the (g-r) color following the prescription given by Bell et al. (2003). 



TABLE 2 

Ultraviolet, Optical, Near-Infrared, Infrared flux densities 



Rad. Ran. 


FUV 


NUV 


u 


g 








J 


H 




3.6 


4.5 


5.8 


8 


24 


70 


160 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


0-75 


0.045 


0.076 


0.19 


0.90 


1.70 


2.33 


2.82 


4.44 


5.28 


4.32 


2.09 


1.32 


1.73 


3.57 


4.02 


>50.29 


>79.72 


75 - 200 


0.008 


0.014 


0.08 


0.43 


0.84 


1.14 


1.35 


1.67 


1.99 


1.54 


0.86 


0.56 


0.60 


0.96 


0.78 


<13.62 


<54.33 


200 - 430 


0.005 


0.010 


0.06 


0.30 


0.55 


0.67 


0.77 








0.45 


0.32 


0.21 


0.29 


0.21 


<4.15 


<22.88 


- 430 


0.058 


0.100 


0.34 


1.65 


3.14 


4.20 


5.01 


>6.18 


>7.37 


>5.94 


3.44 


2.22 


2.56 


4.85 


5.02 


69.85 


167.34 



TABLE 3 

Best-fit parameters and derived quantities of the SED GRASIL models 



Rad. Ran. 

(") 


n 
(Gyr) 


(1O 1O M ) 


Vsch 

(Gyr- 1 ) 


M mo! /M 9 as 


tesc 

(Gyr) 


(10 6 M q ) 


SFR 

(M yr- 1 ) 


TMC 


M 9as 
(10 8 M Q ) 


(1O 1O M ) 


0-75 
75 - 200 
200 - 430 
- 430 


3.7±0.9 
2.8±0.9 
2.2±0.6 
3.9±1.1 


3.1±0.5 
1.4±0.3 
1.0±0.1 
4.8±0.8 


2.4±0.8 
2.1±1.0 
2.3±0.8 
2.3±0.8 


0.49±0.25 
0.52±0.24 
0.71±0.14 
0.49±0.25 


0.030±0.017 
0.051±0.023 
0.032±0.009 
0.031±0.019 


1.3±0.7 
6.9±2.9 
8.9±1.3 
4.4±3.4 


0.64±0.40 
0.19±0.18 
0.07±0.10 
1.07±0.62 


60±14 
346±98 
486±12 
205T121 


2.63±0.62 
0.84±0.31 
0.28±0.26 
4.50±1.00 


2.75±0.45 
1.27±0.30 
0.86±0.62 
4.22±0.63 



Panchromatic view of M94 (NGC4736) 
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Fig. 5.— Left panel: SFR profile of M94. Right panel. sSFR profile of M94. 
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Fig. 6. — M94 (NGC4736) galaxy 0.15-160/xm SEDs of the total galaxy, its bulge, inner disk and outer disk regions. The solid points 
are the observed data and the line is the best fit model. 




R (kpc) 

Fig. 7. — Fit to the surface density profile (shown in Fig. 4) using different analytical components. The oval distortion (inner disk) is 
modeled with a Ferrers bar with n=l. The bulge, outer disk and halo are fit using Miyamoto-Nagai models. 



Panchromatic view of M94 (NGC4736) 15 
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Fig. 8. — Angular frequency of the proposed model. The solid horizontal line is the pattern speed chosen in this paper. 




X (kpc) 

Fig. 9. — Left panel: Snapshot of the SPH simulation showing the spiral arms as a response to the inner oval distortion. On the 
right panel, the solid line shows the gas surface density profile of our model. The different peaks represent the response of the gas to the 
dynamical effect of the oval distortion. The dashed line shows the FUV surface brightness (i.e. a good proxy to the recent star formation) 
scaled to fit in the plot. 



